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Abstract 

The vector and axial pion-photon transition distribution amplitudes are analyzed 
in the Spectral Quark Model. We proceed by the evaluation of double distributions 
through the use of a manifestly covariant calculation based on the a representation 
of propagators. As a result polynomiality is incorporated and calculations become 
rather simple. Explicit formulas, holding at the low-energy quark-model scale, are 
obtained. The corresponding form factors for the anomalous decay tt° — > 77* and 
the radiative pion decays tt^ — ► I vyy are also evaluated and confronted with the 
data. 
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Recently Pire and Szymanowski [T1I2] proposed generalizations of the parton 
distributions to the case where the initial and final states correspond to dif- 
ferent particles. Such objects are termed transition distribution amplitudes 
(TDAs) and are relevant in the analysis of the virtual Compton scattering 
and other exclusive processes. For a recent review of a related topic of gener- 
alized parton distributions (GPDs) see, e.g., [3.4, 5|6|f7] and references therein. 
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Of particular importance are the pion-photon vector and axial leading-twist 
TDAs, V and A, which are the subject of this letter. A quark-model analysis 
of these objects has been undertaken by Tiburzi [8], where the relevant dou- 
ble distributions have been computed. Here we follow similar steps, however, 
instead of using parameterizations [9], we carry out an explicit analytical cal- 
culation all the way down using the Spectral Quark Model (SQM) [10]. This 
model possesses a regularization which allows for a uniform treatment of both 
anomalous and non-anomalous processes, essential in the study. It also en- 
codes the vector-meson dominance and as a result leads to very successful 
phenomenology of numerous processes with pions, photons, or p-mesons. Be- 
low we obtain simple analytic expressions for V and A and the corresponding 
form factors, which hold at the low quark-model energy scale. Variants of chiral 
quark models have been used to obtain information on the pion structure func- 
tion [TTlfTS"] . pion distribution amplitude [T3"jl4lll5|ll6|ll7| . generalized parton 
distributions of the pion [T8"fl9ll2Uf21f22| . and the photon and p distribution 
amplitudes [23J. 

In this paper the considered TDAs are defined as [lj 

-^-e^ + -( 7 (p' )£ )iv;(-^ 7 ^v(^)k 6 (p)}, + =o 

Ztt l l l z t =o 

i p 

=^j^ a ^:p^v ab ^c,t), (i) 

J Z7T 2 2 2 z T =0 

~(e*-q)p^A ab (x,(,t) + ..., (2) 
P + f 

yab = S aby i=Q + ^Sy^ jab = gab ^ + ^ A I=1 , (3) 

where ip represents the iso-doublet quark field, z is a light-cone coordinat^] 
(z + = 0, zt = 0), p denotes the momentum of the pion, / is the pion decay 
constant with / = 86 MeV in the chiral limit, and the photon carries momen- 
tum p' = p + q and has polarization e. Finally, in the asymmetric notation of 
GPDs we use ( = q + /p + and t = q 2 . We consider isovector quark bilinears. 
The isospin decomposition follows from the fact that the photon couples to 
the quark through a combination of isoscalar and isovector coupling, i.e. the 
quark charge is Q = 1/(2N C ) + r 3 /2. The presence of the gauge link operators 
[— z/2,z/2] is understood in Eq. (|TJ2|) in order to guarantee gauge invariance 
of bilocal operators. For brevity, in Eq. (j2J) only the piece proportional to 
(e ■ q)p^ is retained. The part proportional to (e ■ q)q^, indicated by ellipses, 
corresponds to the pion pole term in the t-channel and is not relevant for the 
evaluation of the axial TDA (see [S] for a detailed discussion of the tensor 

1 We use the convention = z° ± z 3 . 
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Fig. 1. The direct (a) and crossed (b) Feynman diagrams for the quark- model eval- 
uation of the pion-photon TDAs. 

structure including the pion pole term). 

The quark-model evaluation of Eq. ( fTll2"l) amounts to the calculation of the 
diagrams of Fig. [U where the + component of the momentum of the hit quark 
is constrained to the value k + = xP + . We denote u as the constituent quark 
mass and use jscu/f as the pion-quark coupling vertex as required by the 
Goldberger-Treiman relation at the quark level. We switch to the Euclidean 
momenta from now o The 1 = and 1=1 components are the following 
combinations of the direct and crossed diagrams of Fig. [TJ V/ =0 = V + V, 
Vj = i = N C (V — V), where the bar indicates the crossed diagram, and similarly 
for A. Since the crossed contributions are related to the direct terms via a 
simple kinematic transformation, we first analyze the direct diagrams. For the 
vector TDA a simple algebra, involving the evaluation of the trace factor, 
yields 



V(xCt) = 2f— S(k-n-x)u 2 

1 ' J (2vr) 4 [k 2 + u 2 }[{k + q) 2 + u 2 }[{k - p) 2 + w 2 ] ' 1 } 

where n 2 = 0, p ■ n = 1, and q ■ n = —(. Our basic technique makes use of 
the a representation for the scalar propagators, which allows for a covariant 
treatment of the 5 function. Thus, 



V(x,C,t) = 2^ 1-0^/^^ (5) 



x da d(3 ^e-^+^-m^?^)-^-^^)^ 



Next, we introduce the variables s = a + j3 + y = j3/s, and z = j/s, shift 
the momentum to k = k! + zp — yq + iXn/ (2s), and then perform the Gaussian 



This is not a limitation. The calculation can be entirely carried out in the 
Minkowski space if the standard boundary condition uj 2 — ► uj 2 — i0 + is incorpo- 
rated. 
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integral over k. Note that since a,/3, 7 > 0, we get < y, z < 1 and also 
y + z < 1. As a result, 



1 1 



a; 



cfy dz6(l-y- z) / ds<5 (x-z- j/C) e - s[uj2+y{l - y)q2+z{l - z)p2+2yzp - q] 





1 1 



h 2 1 dy l dzS ( x ~ z ~ y Q D ( z > y ^ 



(6) 







where 



D{z,y) 



d(l-y-z)uj 2 



lu 2 + y(l — y)q 2 + z(l — z)p 2 + 2yzp • q 



(7) 



is the double distribution for the vector TDA. With our Euclidean vectors 
p 2 = —m 2 , q 2 = —t, and 2p ■ q = m 2 + t for the real photon. Equation (J7J) 
agrees with the result of Ref. [8]. Note that our use of the a-representation 
has led to a direct evaluation of TDAs from the Feynman diagrams. This is a 
very simple alternative to the more involved expansion in moments, typically 
done in similar studies. 



For TDAs the variable ( cannot be assumed to be positive, as the crossing 
symmetry does not relate initial and final states, unlike for GPDs. In general 
we have — 1 < £ < 1. Next, we perform the z integration, which sets z = x—yQ. 
For the case Q > this gives 



V(x,(,t) 



x_ 1— x\ 

6[x(C-x)]f+6[(x-Q(l-x)] J 

) 



dyD(x-y(,y), {( > 0),(8) 



with the first term having the support x G [0,C], and the second x G [CI]- 
For the case ( < we obtain 



/ 



V(x,C,t) 



8tt 2 



l-x 

1-C 



i-c\ 



6[x{(-x)] +9[x(l-x)} 



\ 



dyD{x-y(,y), {( < 0), (9) 



with the support x G [(, 0) for the first term and x G [0, 1] for the second term. 
Thus the support is correct. As expected, the function V is continuous in the 
x variable, with the derivative dV/dx discontinuous at the points x = 0, 

The contribution from the crossed diagram is formally obtained from the direct 
diagram with the replacement andp — > —p—q- Replacing correspondingly x — > 
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£— x and performing the Miinchen transformation: z — > z, y — > 1 — y — z, yields 
the result ^(x, C, t) = V((—x, (, t). The support of the crossed diagram reflects 
the support of the direct diagram. For the case ( > it is x G [—1 + CC]> 
while for ( < we have x G [— 1 + (, 0]. 

The evaluation of the axial TDA proceeds analogously, with the trace yielding 
an additional proportionality factor (2k + q)-e as compared to the vector TDA 
in Eq. (j3J). The result is 



A(x, C, t) = — y ' dy J 'dz5(x-z- y() (2y + 2z - l)D{z, y). (10) 
o o 

In the present case A(x, (,t) = — A(( — x,(,t), with the minus coming from 
the trace factor. 

Expressing the TDAs through the double distributions has the known ad- 
vantage of an automatic verification of the polynomiality conditions which 
ultimately correspond to proper implementation of the Lorentz invariance. 
Polynomiality is manifest from this form, as the moments J dxV(x, (, t)x n in- 
volve integrals with the factor (z+y() n and result in a polynomial in ( of order 
n. We end these general considerations by writing down the double distribu- 
tions in the so called symmetric variables, a, and [3 (not to be confused with 
the previously introduced Feynman parameters), defined asy = (1 + a — /3)/2, 
z = (3. For the vector case we find 



n(<y s) = ^-H-wy m , 

1 ,P) Lui + [l-(a-{3y}qy4 + (3(l-{3)p 2 + (l + a-{3){3p-q' 1 ' 



while for the axial case we find (a + (3)D(a, (3), in full agreement with Ref. 



The theta functions provide the standard integration limits Jq d(3 JjL da. 



The general expressions presented above are formal, as quark models require 
regularization. The choice of a finite regularization is far from trivial, as re- 
quirements of the Lorentz and gauge invariance, as well as preservation of chi- 
ral Ward identities and in particular anomalies (crucial in the present study) 
impose severe and tight constraints. An elegant way of imposing a regulariza- 
tion which obeys these requirements is achieved in the Spectral Quark Model 
(SQM) |l0j . In this model, developed in the spirit of the early work of Efimov 
and Ivanov [24J, the quark mass u is treated as a spectral parameter of a 
generalized Lehmann representation, which is integrated along a suitably cho- 
sen complex contour C. The chirally symmetric effective action constructed in 
Ref. HOI reads 
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V[U,v,a] = -iNcj dujp(uj)Ti\og (ifi - ujU Fj - {j + ^75)) , (12) 
c 

where the trace for a bilocal (Dirac- and flavor-matrix valued) operator A(x, x') 
is given by TyA = J d i xti(A{x ) x)) with tr denoting the Dirac trace and () 
the flavor trace. The matrix U 5 = e i7sT,7r / / = \{l + 75)^ + ±(1 - 75) 17+, while 
U = e in ' T 'f is the flavor matrix representing the pseudoscalar pions in the 
nonlinear representation. The symbols v M and a M represent external vector 
and axial currents. In Eq. (fT2l the spectral density p(u) acts eulator. 
Actually, the expressions for one-quark-loop observables in SQM are obtained 
from the preceding expressions by integrating over uj with the spectral density 
p{oS). For a generic observable F we have 

F SQM = J dup{u)F{uo), (13) 

c 



where F(u) is the quark- model result with the quark mass set to w. In the 
meson-dominance version of SQM [10J we have 



p(u) = pv(u) +ps{u), (14) 

1 1 1 1 48n 2 (qq) 

PvM = —--. , PS{V) 



2Txiuj(l-Auj 2 /M^fl 2JHby 1 ' 27rziV c M|(l-4cu 2 /M|) 5 /2' 

where Ms = My = m p is the p- meson mass [25J . The contour C for the spectral 
integration over u and other details are given in Ref pi)]. SQM generates the 
monopole pion electromagnetic form factor [10]. Interestingly, the model has 
the feature of the analytic quark confinement, i.e. the quark propagator has no 
poles, only cuts, in the complex momentum plane. Moreover, the evaluation 
of low-energy matrix elements in SQM is very simple and leads to numerous 
results reported in Ref. [TU], in particular for the pion light-cone wave function, 
the pion structure function, the generalized parton distributions of the pion 
[2T] , low energy chiral Lagrangeans [25j and the photon, and p-meson structure 
functions [231. 



Let us assume for simplicity from now on the chiral limit, p 2 = 0, and the real 
photon, (p + q) 2 = 0. The evaluation of the spectral integral is straightforward 
using the results of Ref. [25J. Then the double distribution for the vector TDA 
in SQM assumes the simple form 



1 — y — zj 

D SQM (z,y) = I dup(u)D(z,y) = (15) 



c (1 ^ 



4y(l-y-z)t 
M. 



1 v 



Remarkably, completely analytic expressions for the TDAs follow, fulfilling all 
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a priori Lorentz, chiral and gauge invariance constraints. To our knowledge 
this is the first explicit calculation of TDAs in a regularized chiral quark 
model, and we list the results in some detail as their form may guide the used 
parameterizations of TDAs. We only show the results for the case ( > 0, as 
for C < they have a similar character. We find for the vector TDA 



V SQ m(x, C, t) = ^ (0[x{C - x)] (xi + y ) + ~x)(x- ()) X 2 ) , (16) 
where 



= 2(x-l)[3(C-l)M^ + t(x-l) 2 ] 
[(C-l)M 2 +t(x-l) 2 ] 2 ' 
(x(C - 2) + Q (3M 2 (C ~ 1)C 2 + t ((C 2 + 8C - 8) x 2 + 2(4 - 5QCs + C 2 )) 

((C-i)M 2 +t(*-i) 2 ) 2 (c 2 + ^^j 

Some special values are 

247T 2 (My + t(C - 1)) 

V SQM (x, C, 0) = JL La;( C - x)]^ + - x)(x - C)]^ " ^ 



2 • 



V fr t) - ^(3M 2 -t(l-x) 2 )(l-x) 
Wx ' °' f) ~ 24vr 2 (My - t(l - x) 2 ) 2 • (18) 

The functions 87t 2 Vsqm(^, C> are shown at the top of Fig. [2]for two values of 
t and for several values of £, both positive and negative. The integration over 
x produces the (^-independent (as required by polynomiality) form factors, 



M 2 



24tt 



2 



Ml-t t 



j d xVi Q M(x,(;t)= J dx (Vsqm + ^SQM 

V ' J 

J dxVi^(x, C; t) = J dxN c (V S qm - ^ SQ m) = 0. (19) 
The isoscalar form factor is related to the pion-photon transition form factor, 



^.( ( )^/^».J 37 ( 1 + JL + ^...), (20) 
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where the factor of 2 comes from the fact, that either of the photons can 
be isoscalar. We read out the corresponding rms radius to be (r 2 )]/^* = 
\/5/M y = 0.57fm for My = 770MeV. Equivalent^, one may use the slope pa- 
rameter = f t F^ ir {t)/F^ ir {t) SQM gives b n = 5/(6M£) = 1.4 GeV" 2 , 
in very reasonable agreement with the experimental values quoted by the 
PDG [25]: K = (1-79 ± 0.14 ± 14)GeV~ 2 originally reported by the CELLO 
collaboration [27], b n = (1-4 ± 1.3 ± 2.6)GeV" 2 from Ref. [28], or b w = 
(1.4 ± 0.8 ± 1.4)GeV" 2 given in [29]. 

The vector form factor is also related to the Fy form factor appearing in the 
radiative pion decays, tt^ — > l^vtf (for a review see e.g. [30] and references 
therein). With the assumption of CVC it is related to the the isoscalar form 
factor in the following way: 



Fv(t) = ^ dxV J = = ^ 1 + -^ + ...). (21) 



— J dxV =wl 1 + M? v - , 

The premultiplying factors follow the assumed convent iond" 3 !. The value at 
t = (for / = 93 MeV) is F v (0) ~ 0.026, as listed in [26]. The experimental 
data fall one standard deviation below this CVC prediction, with Fy(0) = 
0.017 ±0.008 



For the axial TDA the corresponding expressions are (for ( > 0) 



My 



A S qm(x, C,t) = ^ (9[x(C - x)] X 3 + 0[(l -x)(x- C)]x X2 ) , (22) 
where X2 is given by Eq. <HM and 



_ (C - lfM v + t(x - l)(5z -2)(C- 1)M 2 + t 2 {x - lf{2x - 1) 
X3 ~ t(t(x-l)2 + (C-l)M 2 ) 2 + 

1 

t (M 2 + ^Of 2 ( t (x - l) 2 + (C - 1)M 2 ) 2 
x -(C - 1) 2 ( 3 M V + t(C - 1)C 2 ((C - 6)x 2 + 7(x - 2C) M 2 + 
t 2 (16a; 4 - 24(x + l)(x 3 + G(x(x + 6) + 1)(V - {5x(x{x + 3) - 1) + 1)C 3 ) 

The special values are 



(23) 



3 We follow |26|30j . The structure dependent amplitude is given by Msd = 
ieG cos c u„'fij,(l — ^^)v e e* u MV /^.m^, where the hadronic contribution is given 
by AP" = Fy(t)e^ Paqf3 - iF A {t)[q v (q» + p») - gTq • (q + p)}- 
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-0.2 L w -0.2 L 

Fig. 2. Top: vector TDA 8Tr 2 V S Q M (x, (,t) for t = (left) and t = -0.4 GeV (right) 
plotted as functions of x for several values of Q: —1, —2/3, —1/3, 0, 1/3, 2/3, and 1. 
The value of £ can be inferred from the position of the cusp for ( > or from the 
support for ( < 0. Bottom: the same for the axial TDA 87t 2 Vsqm(:£, C> The vector 
curves are normalized to 1/2, and the axial curves to 1/6. The presented calculation 
is made in the Spectral Quark Model with My = 0.77 GeV in chiral limit and for 
the real photon. 

AsQM(0,C,t)=4sQM(l,C,<) = 0, 

A (ttt ,_ CM v (3M v + t(C-l)) 

^s Q m(c, c, *) - 247r2(M , +t(c _ 1))2 - owe, C, *), 

-4sqmO,C,0) = 



•571 



g[x(C-x)] ^±^ +6)[(l-x)(x-C)]^-^ 



-4sQM(a;,0 



_ My (3M 2 



- X 



(1 — x)x 



24vr 2 (M y 



(i-C) J' 

xV SQM {x,0,t). (24) 



The form factors are 



dx (v4 SQM + A S q M 
dxA SQM( x '(,t)= J dxN c (A SQM -A 



SQM 



MyiV c 



log 1 



(25) 



247T 2 t 

The bottom part of Fig. [2] displays the function 87t 2 Asqm{x, C t)- 

The axial TDA is related to the axial-vector form factor measured in the 
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radiative pion decays 7r =1= — > / ± z/;7 via integration over the x variabld_], 

FM = ^ J dxA .^ c, () = ^ {, + _±, + JL + . . .) , (26) 

where the premultiplying factors are the same as in Eq. (1211) . The value at 
t = is F A (0) = F v (0) ~ 0.026, which is a factor of 2 larger than the 
experimental number Fa{0) = 0.0115 ± 0.0005 [26J. The same conclusions 
were reached in Ref. [8j. The predicted ratio Fa(0) / F v (0) = 1 compares to the 
experimental number of Fa(0)/Fv(0) = 0.71^2 within one standard deviation. 
The t-dependence of the form factors is presented in Fig. [31 We note long tails 
due to a log(— t)/t behavior at large t. The axial form factor (dashed line) lies 
above the vector form factor (solid line), which in turn lies above the monopole 
vector- meson dominance form factor, drawn as reference (dotted line). 

It is interesting to analyze our results in the light of Chiral Perturbation 
Theory in the large-iV c limit [31], where F A (0) = 4(L 10 + L 9 )\/2m^/ = 
0.012 ± 0.008 for L 9 = (6.9 ± 0.7) • lO" 3 and L w = (-5.5 ± 0.7) • 10~ 3 (in 
fact the values of the low-energy constants L 9 and L 1Q are determined from 
the pion electromagnetic form factor and the radiative pion decays). In the 
large- N c limit one imposes QCD-motivated short- distance constraints regard- 
ing the asymptotic falloff both for the difference of two-point vector minus 
axial correlators (the Weinberg sum rules) as well as the axial and electro- 
magnetic form factors at large t (see, e.g., [32]). In the single resonance ap- 
proximation (SRA) one gets the same value for Fy as above and since Lg RA = 
-4LSR A /3 = p/2M v a ratio F A (0)/F v (0) ~ 1/2 = 0.35 which produces 
a significantly lower value for the axial form factor, Ff RA (0) ~ 0.01. In SQM 
one has L S 9 QM = -2L S ^ M = N c /(A8n 2 ) [25] and hence F^ QM (0) = 0.026 in 
agreement with Eq. (126]) . The mismatch in the L\q/ Lg ratio in SQM and SRA 
stems from the absence of an explicit axial meson contribution in SQM which 
induces the violation of the second Weinberg sum rule pointed out previ- 
ously [I0] and generating a Fa(0) about twice the experimental number. This 
feature is common to all known local quark models. The SQM axial radius 
and the large- iV c SRA result coincide, {t 2 )a = 3/My, although at large t SRA 
yields ~ 1A which is slightly more convergent than our result (125]) . 

which contains an additional log(— t). It would be interesting to pursue the 

4 This yields, according to Eq. ([2]), the quark contribution to the axial current from 
the 7^75 vertex only, which by itself does not convey the chiral Ward-Takahashi 
identity [10] at the quark level. The additional pion pole contribution — 2uq t1 ^ /q 2 
must be included in the vertex. Although this is essential to preserve the identity 
("f(p' , e)\J^~ \tt + (p)} = ie(e* ■ q)V2f n , and therefore for an unambiguous identifi- 
cation of Fa{P) from the matrix element, for the tensor structure in Eq. (f2]) such a 
contribution turns out to cancel. Hence in the evaluation of the axial TDA we may 
retain only the pieces proportional to in Eq. ([2]), obtained with the 7^75 vertex. 
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Fig. 3. Form factors i^^ 77 * (i)/.F 7r _ v77 * (0) (solid line) and -Fa(*)/-^a(0) (dashed 
line). The dotted line shows as reference the electromagnetic form factor which in 
SQM is a monopole F% m {t) = M^/(M^ - t). 

present calculation to the nonlocal chiral quark models where both Weinberg 
sum rules are known to hold [33]. Finally, the large t-behavior is subjected to 
the QCD logarithmic radiative corrections. 

Actually, the results obtained above hold at a low energy scale of the quark 
model. In Ref. [H] an estimate of this scale based on the momentum sum 
rule has been given. For the present model or other local models, such as 
the Nambu-Jona-Lasinio model, the scale turns out to be very low, around 
320 MeV. An independent estimate based on the pion light-cone distribution 
amplitude results in a very similar estimate [17]. For that reason, evolution is 
crucial for the case of DAs of PDFs, and undoubtedly will also be important for 
the present case of TDAs as well as the vector and axial form factors at large 
momenta. The results obtained so far simply provide the initial conditions 
for the QCD evolution, which at the leading order can be made with the 
usual ERBL equations [3lf35] . A detailed study of the evolution issues will be 
presented elsewhere. 

In conclusion, we have presented an explicit quark-model study of the vector 
and axial leading-twist pion-photon TDAs. The results are analytic, which 
allows us to gain insight into the possible forms of non-perturbatively gener- 
ated TDAs. Such predictions for TDAs or GPDs are scarce, and frequently 
one only makes guesses subject to the constraints coming from form factors, 
polynomiality, etc. Our method conforms to the Lorentz and gauge invariance, 
preserves the chiral Ward-Takahashi identities as well as satisfies anomalies, 
crucial in the study of the VAA processes. Since we proceed via the double dis- 
tributions, our TDAs automatically satisfy the polynomiality constraints. The 
used technique of calculation, which takes advantage of the a-representation 
of propagators and thus is manifestly covariant, makes a direct use of the 
Feynman diagrams and produces the results in a few straightforward steps. 
No expansion in terms of moments is necessary. This technique is applica- 
ble also in other calculations of this kind: PDFs, GPDs, as well as for other 
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models, including the non-local models, where the quark mass depends on the 
virtuality. Our results correspond to a low-energy quark-model scale. After 
suitable QCD evolution the obtained results may be used in the studies of 
the virtual Compton scattering and other exclusive processes involving pions, 
photons, and the weak gauge bosons. 

We thank Brian Tiburzi for helpful e-mail exchanges and for pointing out a 
mistake in one of our formulas. 
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